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I , Abstract. Fractional, anomalous diffusion in space-periodic potentials is investi- 

"t^ ' gated. The analytical solution for the effective, fractional diffusion coefficient in an 

"^ ■ arbitrary periodic potential is obtained in closed form in terms of two quadratures. 

1 t I This theoretical result is corroborated by numerical simulations for different shapes of 
^ • the periodic potential. Normal and fractional spreading processes are contrasted via 
Vh ■ their time evolution of the corresponding probability densities in state space. While 

1-^ I there are distinct differences occurring at small evolution times, a re-scaling of time 

P^ ■ yields a mutual matching between the long-time behaviors of normal and fractional 

O diffusion. 
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CN . 1. Introduction 

o ■ 

\0 [ In 1905 Pearson proposed what we now know as a random walk pTj. For a one- 

^ I dimensional system the problem can be formulated in the following way: A particle 

c^ I jumps at each point of time from its current position x to the position x + Ax with 

G I probability p, or x — Ax with probability 1—p. The approach towards diffusion theory, 

Ti \ pioneered by Einstein, relies on postulates very similar to the ones for the random 

Q I walk, leading to the same results when the jump width Ax — > ™. However, in 

J-^ I many situations, the assumptions used by Einstein and Pearson do not hold; one of 

.^ I such examples is the transport of charge carriers in amorphous semiconductors when 

rS I exposed to an electric field. 

c^ I Sixty years after Pearson, in 1965, Montroll and Weiss introduced the theory 

of continuous time random walks (CTRW) jH]. It was applied to the transport in 
semiconductors in works by Scher and Lax [^, and Scher and Montroll jS]. Due to its 
historical importance and vivid clarity we recall here the definition of the CTRW as given 
by Scher and Montroll [5]: In our model we postulate our material to be divided into a 
regular lattice of equivalent cells, with each cell containing many randomly distributed 
localized sites available for hopping carriers. Carrier transport is a succession of carrier 
hops from one localized site to another and finally from one cell to another. We define 
the hopping time to be the time interval between the moment of arrival of a carrier into 
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one cell and the moment of arrival into the next cell into which it lands. The random 
distribution of sites and hence the disorder of an amorphous material is incorporated 
into a hopping-time distribution function iP{t). The appropriate distribution iP{t), 
leading to the agreement with the experiments, was shown to possess the power-law 
form: ^(r) oc r~^~" with a G (0, 1) JSHEIIZI- For this range of the fractional exponent 
a all the moments of the distribution 'iP{t) diverge and the corresponding process has 
no characteristic time scale, thus exhibiting the phenomenon of aging. As a result, the 
process undergoes subdiffusion [31011101, i-e., the mean square displacement grows in 
the absence of an external force slower than linearly in time, {6r'^(t)) ~ t" (0 < a < 1). 
In the original study of the fractional transport in the context of anomalous 
transport in semiconductors an ensemble of carriers executing a random walk, when 
biased by an electric field, was studied [3 El- In the present paper our focus is 
different: We instead address the problem of the carriers executing the random walk 
in a spatially varying periodic potential. This situation is representative for many 
applications occurring in areas such as in condensed matter physics, chemical physics, 
nanotechnology, and molecular biology, to name but a few. For those applications it 
is of utmost importance to account for the spatial variation of the transport process 

[iiiiiaiisiiiiiiiaiini. 

Our work is is set up as follows: In section |21 we propose the model and define 
the theoretical and numerical problem. In section JHl we recall some prior results about 
the biased CTRW [T71 and a CTRW proceeding in a washboard potential [HI. In 
section |3] the formula for the effective fractional diffusion coefficient in a periodic 
potential Uq{x) = Uo{x + L) with period L is derived and the theoretical result is 
corroborated by numerical simulations of the CTRW for different shapes of periodic 
potentials. Finally, we address the problem of particles spreading anomalously in a 
periodic potential also in the light of the time evolution of the space probability density, 
as compared to the case with normal diffusion. 

2. Set up of the model 

Following the general picture of the CTRW we introduce a one-dimensional lattice 
{xi = iAx} with a lattice period Aa; and i = 0, ±1, ±2, . . . After a random residence 
time r a particle at site i hops to site i±l with a probability q^ (see figure ^a)): The 
sites here correspond to the cells in reference |3] . The random time r is extracted from a 
residence time distribution iP{t). A suitable possible choice for iP{t) is a Mittag-Lefiler 
distribution defined by 

^.(r) = -;fi?.(-(^.rr), with E^{-{u^,rr) = f^t^^^^ . (1) 

dr ^-^ i (na + 1 

n=0 ^ ^ 

The quantity z/j is the time-scaling parameter at lattice site i. 
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The CTRW with the Mittag-Leffler residence time density can be described through 
a fractional master equation for the site populations Pi{t) fSlQSl; i-^., 

D':Pi{t) = /,_i P,;_i(t) + (7,+i P,+i(t) - (/, + Qi) P,{t) , (2) 

with the Caputo fractional derivative D'^ 30^ on the left-hand side defined by 

1 /■*,., 1 d 
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(3) 
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r(l - a) Jo "if- t'Y dt' 
The quantities fi = qfv'^ and Qi = q^i'" in the fractional master equation © are referred 
to as the fractional forward and backward rates. Using the normalization condition for 
the splitting probabilities, i.e., qf + g" = 1, one obtains that 

^^ - U/ifi + 9i) , C = 9r/{fi + 9d , (4) 

{f^ + 9^Y^''■ (5) 

For an arbitrary shaped potential landscape U{x) the fractional rates can be chosen as 

/, = (/€ JAx^) exp[-/?(f/,+i - U,)/2] , (6) 

g, = (kJAx^) exp[-/3(f/,_i - f/,)/2] . (7) 

Here Ui = U{iAx) and f3 = l/ksT is the inverse temperature; Kq is the fractional 
free diffusion coefficient with dimension cm^s""". The form ©-© of the fractional 
rates ensures that the Boltzmann detailed balance relation is satisfied, i.e., fi-i/9i = 
exp[P{Ui-i — Ui)]. The lattice period Ax must fulfill the condition U"{x)Ax <^ 2U'{x) 
in order to ensure the smoothness of the potential. In the case of a periodic potential 
this implies in particular that the lattice step size Ax is much smaller than the potential 
period L, Ax <^ L. Furthermore, in order to recover the continuous limit addressed 
below, the condition \P{Ui±i — f/j)| <^ 1 must be obeyed. 

In the space-continuous limit the CTRW with the Mittag-Leffler residence time 
density can be described through the fractional Fokker-Planck equation [HI CHI EDI EI] , 

d U'(x) d^ 



D:P{x,t) 



dx 



Vo 



+ Kr 



dx"^ 



P{x,t). 



(8) 
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Figure 1. (color online) CTRW in a periodic cosine-potential. Two different 
possibilities to introduce the one-dimensional lattice: (a) The lattice period Ace is 
much smaller than the potential period L. The particle at site i hops to site i -I- 1, or 
i — 1, respectively. (b)The lattice step size Ax is equal to the potential period L and 
the lattice sites are centered at the potential minima. The particle performs hops from 
one potential minima to one of the neighboring ones. 
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Here, P(x, t) is the probability density and a prime stands for the derivative with respect 
to the space coordinate. The quantity rj^ denotes the generahzed friction coefficient, 
possessing the dimension kgs"~^. It is related to the bare fractional anomalous diffusion 
coefficient /€„ through ?7q,Kq, = k-oT. 

With this material at hand we have defined our theoretical problem as well the 
numerical procedure. In the simulations of the CTRW we use for a G (0, 0.8] a Pareto 
residence time distribution, i.e., 

^.(r) = -^P.(..r), with ^«(-^-) = ir^^T^^yiA;^ > (9) 

instead of the Mittag-Leffier one, as for every < a < 1 the long time behavior of 
the system is determined solely by the tail of the residence time distribution J22I- For 
a > 0.8 the Mittag-Leffier density (0) is employed. For a = 1 the latter one transforms 
into the exponential distribution, covering the regime of normal overdamped Brownian 
motion. The spatial lattice step in our simulations is Aa; = 0.001, measured in units of 
the spatial period L. The energy is measured in units of the potential amplitude A, and 
the time unit is set as tq = (rjaL^/Ay^". For a detailed description of the algorithm 
for the numerical simulations and of the employment of the Pareto or Mittag-Leffier 
distribution, we refer the readers to the comprehensive work in reference |2Hj . 

3. Biased CTRW and CTRW in a washboard potential 

3.1. Biased CTRW 

The anomalous diffusion that is biased by a constant external force F is a well established 
phenomenon found in many different systems. For the biased CTRW the fractional rates 
(0)-© become site-independent, fi = f and gi = g, as Ui±i — Ui = iFAx. From the 
fractional master equation Q one finds then the solutions for the mean particle position 
and for the mean square displacement [^, 

(x(t)) =(x(0)) + ^^t^r, (10) 

V[a + 1) 

{5x\t)) = {6x\0)) + ^^il±3lt» 



+ 



r(a + l) 

1 



X Ax'if-gyt"'. (11) 

The solutions of the corresponding fractional Fokker-Planck equation are in the 
same form of the ones for the fractional master equation; i.e., 

(x(t)) =MO)) + -^.T^, (12) 

r]a T{a + 1) 

{6x\t)) = {6x\0)) + 2k ^" 



+ — 



r(a + l) 

1 



r(2a + l) r2(a + l) 



t'" . (13) 
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The comparison of the solutions (fTI]|) to (jT^ and (fTTj) to (jT^ gives, 

Ax(/ - (7) = F/r^, and Aa;2(/ + (?)/2 = /., . (14) 

The latter equations define the anomalous current and the anomalous diffusion 
coefficient through the fractional rates / and g, and are of the same form as the 
corresponding relations for the normal Brownian diffusion, determined through the 
corresponding escape rates. 

If at a given temperature T the system is close to thermal equilibrium, the mean 
square displacement in the absence of an external force and the average displacement 
induced by a bias F 7^ are related through the generalized Einstein relation jHl 1211 123, 

(5^'W)U=o = ;^ (a:(t)-x(0))|^. (15) 

Note that equation (jl5|) is strictly valid only in the linear response regime, which is 
approached when when F — i> 0. It then leads to the generalized ffuctuation-dissipation 
theorem 

t^a = iPVa)-' . (16) 

3.2. CTRW in a washboard potential 

Solving in the stationary limit the fractional Fokker-Planck equation (jH} for a biased 
periodic potential U{x) = Uq{x) — Fx, one finds for the mean particle position |18j, 

{x{t)) = {xm + -^l^^t-. (17) 

T{a + 1) 
The anomalous current Va{F) in the washboard potential is then given by a generalized 
Stratonovich formula, put forward in reference JH], i-e., 

K^L [1 - exp{-f3FL)] 



VaiF) 



/,^ dx /;+'' dy exp(-/3[f/(x) - Uiy)]) 



/o 

In analogy to the case with normal Brownian motion, in order to study the fractional 
diffusion in a periodic or washboard potential, it would seem natural to choose the lattice 
period Ax to be equal to the space period L and the sites to be centered at minima, 
as illustrated in figure ^b). In this case the fractional rates, that we mark for such a 
lattice with fj and gj, are independent of the site j, fj = f and gj = g. It was proved 
in reference [18j that considering the CTRW in the lattice {xj = jL}, the asymptotic 
solution (t -^ 00) for the mean square displacement in a tilted periodic potential is 
of the form as the solution (fTT|) for the biased CTRW, as the equation (fTTj) is of the 
same form as equation (fTUj) . The fractional rates / and g, however, are no longer given 
by equations (jEI)-© and the model does not provide their explicit dependence on the 
potential. 

Whereas in the washboard potential Ax(/ — g) is equal to the generalized 
Stratonovich current Va, one could expect that Ax^(/ + g)/2 follows a generalized 
formula for the effective diffusion coefficient in a tilted periodic potential [23 I2Z], in 
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correspondence to equations (jT^ because the problem can be mapped onto the case 
with a constant bias. However, in the long time limit the ballistic term oc t^" prevails 
over the term proportional to t" and the effect of the latter one is negligible: The ratio 
between the mean square displacement and squared average coordinate depends in the 
asymptotic limit only on the fractional exponent a and obeys the same result as for the 
biased CTRW |H1QH1- The term proportional to t" becomes relevant for t -^ oo only in 
the limit a — > 1, leading to the normal diffusive behavior, or for F— ^Oas/ — ^^0, 
i.e., for a periodic potential. 

4. Fractional diffusion in a periodic potential 



4.1. The mean square displacement 

In this section we present our results for fractional diffusion in spatially varying, periodic 
potentials. We start from the expression of the mean square displacement for the particle 
in the periodic potential Uq{x). For zero tilting the fractional rates / and g become 
equal and the ballistic term occurring in equation (fTT|) . reformulated for a washboard 
potential, thus disappears. Therefore, the asymptotic mean square displacement now 
reads, 

{5x\t)) = {6x\0)) + Ax\f + g) — ^ . (19) 

i (a + 1) 

This equation is confirmed by the numerical results, depicted in figure El for various 

values of the fractional exponent a, which is equal to the slope of the numerically 

evaluated curves on the logarithmic scale. 

Correspondingly, the asymptotic solution of the fractional Fokker-Planck equation 

for the periodic potential Uq{x) can be written in the form analogous to the case of the 

fractional diffusion in the absence of force, i.e., 

fr^TTT ' (20) 

i a + 1 
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Figure 2. (color online) The mean square displacement computed numerically for 
the CTRW in the cosine-potential U{)[x) — Acos(27rx/L) for different values of the 
fractional exponent a. The re-scaled temperature is set at k^T j A = 0.5. 
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This equation defines the effective fractional diffusion coefficient k. 



(eflF) 
■a 1 



4-"i = r(. + i)mnM^l)H^£!M. (21, 

Also the exphcit expression for the fractional rates to hop from one minimum to one of 
the neighboring minima follows from equations (fT^ and (j^UI) . as for f = g 



.f ) = MZ±1) ^ Ax^f ^ Ax^g 



^(en) ^ ^_^ ^ ^^,f ^ ^^.^ _ (22) 

4-2. Effective fractional diffusion coefficient 

Next, a useful analytical expression for Ka in a periodic potential can be derived from 
a generalized Einstein relation. Equation (fTH|) is valid in the linear response regime also 
for any periodic potential Uq{x), as the problem of fractional diffusion in a periodic 
potential can be mapped onto the force free case |^. In doing so we can write 



^(eff) - 1 hm !^^ - 1^!^^ 

^- -/5f^o F ~ (3 dF 



(23) 



F=0 



where Va{F) is given by the generalized Stratonovich formula derived in (fTHj) . As a 
central result we thus obtain the following closed, exact analytical expression for the 
effective fractional diffusion coefficient, reading 

K^f^ = J ^ . (24) 

L 2 J^ dx exp [(3Uo{x)] J^ dy exp [-f3Uo{y)] 

This expression is valid for an arbitrary shaped, unbiased periodic potential Uo{x). It 
reduces for a = 1 to the corresponding formula for the normal diffusion in a periodic 
potential, first derived by Lifson and Jackson in reference [21] and independently again 
in references \^'2\ ES] ■ Our new result therefore provides the generalization for fractional 
diffusion processes which are anomalous. 

The behavior of equation (J2l|l versus re-scaled temperature k-sT/A is illustrated in 
figure El^a) for the following periodic potentials, depicted with figure Efb): 
(i) a cosine potential, 

[/g(i) (x) = A cos{2nx/L) , (25) 

(ii) a double hump potential, 

t/^2) (^) ^ ^^^ [cos(27rx/L) + cos(47rx/L)] , (26) 

with the coefficient ai = 16/25, and 
(iii) a ratchet potential, 

f/g(3) (a;) = A[a2 sm{2TTx/L) + a-s sin(47rx/L)] , (27) 

with 02 = 85/(21v^), as = 25/(21^/21)- The coefficients Oi, 02, a^ are chosen such 
that the potentials (j^^jl - ljTTj) have the same amplitude A. The theoretical curves are 
confirmed by numerical results, depicted in figure El with symbols. The anomalous 
diffusion coefficient is computed as defined by equation (J2H). As the ratio Ka /na < 
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1, one can conclude that, analogously to the normal case, the effect of any one- 
dimensional non-biased periodic field is to suppress the macroscopic anomalous diffusion 
coefficient compared to the value in the absence of force [STj . A possible enhancement 
may be expected in presence of time-dependent, periodic landscape modulations as 
demonstrated for normal diffusion in references [S2 ESI IHE]- Furthermore, it is to 
be noticed that the ratio Ka / i^a does not depend on the fractional exponent a and 
moreover, the shape of the periodic potential Uq{x) has only a small influence, as one 
can see by comparing the theoretical curves in figure El (a) (note also references IHZIEHI)- 



4-3. Probability density: Anomalous versus normal 

In the previous section it was demonstrated that the effective fractional diffusion 
coefficient in a periodic potential is of the same form as the the Lifson- Jackson formula 
for normal diffusion. This represents a further element of the formal analogy between 
fractional and normal diffusion, besides e.g. the validity of the generalized Stratonovich 
formula ()18|) ^Hl and the fact that the stationary reduced probability density is the same 
for both cases [SB]- Here, we present additional results which support and corroborate 
this formal analogy further. We notice that in the absence of a bias, all the odd moments 
of the probability density are identically zero both for normal and fractional diffusion. 
As for the second moment, upon introducing the re-scaled time, 
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Figure 3. (color online) (a) Effective anomalous diffusion coefficient Ka in a periodic 
potential versus the re-scaled temperature k^T /A. The quantity Ha is re-scaled 
by the corresponding free fractional diffusion coefficient Ka. The theoretical curves 
obtained from equation (|24|l (lines) are compared to the numerical results (symbols). 
The different periodic potentials used are given by equations (I25II - H27II . For each 
potential and at given temperature the numerical points are computed for some values 
of a within the interval a € [0.1, 0.9]. (b) A comparison among the different periodic 
potentials used for the numerics, see in (P51) - (|77|l : (1) the cosine potential C/q {x); (2) 
the double hump potential U^ (x) — (2ai — 1); (3) the ratchet potential Uq (x). 
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it follows from equation (j2Uj) that the mean square displacement (in units of L^) formally 
coincides with that of the normal diffusion case, [{5x'^{t')) — ((5x^(0))]/L^ = 2T't', 
independently of the fractional exponent a, wherein T' = ksT/A, with A the potential 
amplitude, is the re-scaled temperature. The study of the time evolution of the 
probability density is illustrated with the example in figure ID choosing the times t for the 
anomalous diffusion process and the corresponding times t' for normal diffusion, so that 
they satisfy equation (j^Hj) : The probability densities for anomalous diffusion (continuous 
lines) and normal diffusion (dashed lines) processes are barely distinguishable from each 
other for sufficiently long evolution times. 

In clear contrast, however, appreciable differences between the normal diffusion 
coordinate density P{x,t') and the anomalous coordinate density P{x,t) emerge for 
small times. This is best detectable by comparing the reduced probability density, 
mapped onto a single spatial period, 

P{x,t) = '^P{nL + x,t), neZ, (29) 



as done in figureEl In the normal case (figure El left) the two initial maxima at x = and 
x/L = 1, due to the initial conditions P{x, 0) = S{x), move toward the center and finally 
merge into the asymptotic stationary density (solid line) Pst{x) = Af^^ exp[—i3Uo{x)], 
where A/" = /^ dx' exp[—i3Uo{x')] is a normalization factor and Uq{x) = Acos{27!-x/L). 
On the other hand, in the anomalous case the two initial maxima gradually disappear, 
while a new peak grows at x/L = 0.5 and evolves into the stationary density Pst{x). 

Moreover, as soon as the process is biased by an external finite force, F 7^ 0, a 
qualitative difference arises in the time evolutions of the probability densities of the 
anomalous and the normal processes in the long time limit as well, see also Ref. |23j . 
This is true even for small values of F in the linear response regime, as one can defer from 
figure ini All this indicates a profound difference between a fractal diffusion dynamics 
that is based on the fractal Brownian motion introduced by Mandelbrot and van Ness 
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Figure 4. (color online) The time evolutions of the probability densities characterizing 
the anomalous and normal diffusion processes in the periodic cosine potential (|25|l . The 
re-scaled temperature is kBT/A — 0.5 and the fractional exponent is a = 0.5. The 
anomalous probability density P{x, t) cannot be distinguished from that of the normal 
case, P{x,t'), once the time has been re-scaled according to equation (|28|) . Similar 
results are obtained for other values of a G (0, 1) (not depicted). 
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Figure 5. (color online) The different small-time evolutions of the normal (left) and 
anomalous (right) reduced probability densities P{x, t) and P{x, t') defined by equation 
(|^ . in the cosine potential (^5)1 . Curve labels (1), (2), (3) and (4) represent increasing 
values of re-scaled time t' = 0.01, 0.02, 0.04, 0.11 for the normal case and of time t for 
the anomalous case, related to t' through equation l|28|l . The solid line represents the 
theoretical stationary solution. The re-scaled temperature is k-gT/A = 0.5 and a = 0.5 
for the anomalous process, as in figure 01 
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Figure 6. (color online) The time evolutions of the probability densities characterizing 
the normal (above) and anomalous (below) diffusion processes in a tilted cosine 
potential U{x) = A cos{2'kx/L) — Fx. The re-scaled temperature is k-oT /A — 0.5 and 
the fractional exponent is a = 0.5, as in figure 01 The tilting force is i^ = 0.1 x Fcr, 
where Fcr — 2tiA/L is the re-scaled critical bias, corresponding to the disappearance of 
potential minima. For sufficiently small times the probability densities of the normal 
and anomalous processes are very similar. However, at larger times (in the long time 
limit) the maximum of the density for normal diffusion moves with the directed current. 
In contrast, the mean square displacement of an ensemble of particles undergoing 
fractional diffusion is dominated by the ballistic contribution and the typical stretched 
spreading in the direction of bias is observed, while leaving the maximum of the density 
near the origin. 
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PH] and the fractional diffusion based on the CTRW jH]. The time evolution of the 
density of an ensemble of particles undergoing normal diffusion can be interpreted as a 
superposition of a translational motion and a spreading of the initially localized density. 
In this case one observes the global maximum of the probability density moving in the 
direction of the external bias (figure IHl (above)). Instead, in the anomalous case, only a 
spreading of the initial density takes place, resulting in a long tail in the direction of the 
bias. The global maximum of the density remains close, however, to its initial position 
(figure El (below)) [SB]- This intriguing behavior is related to the presence of a ballistic 
contribution proportional to f^^ in the mean square displacement (see equation (fTTj)). 
We remark that for a close to one and for small values of external bias F, at small 
times the term oc t" can prevail the ballistic term. However, in the long-time limit 
the ballistic term takes over and always dominates. The latter remark may be relevant 
for experimental studies. It in addition also provides a crucial test that allows one to 
distinguish between fractal and fractional Brownian motion on a practical level. 

5. Conclusion 

With this work we investigated anomalous diffusion whose dynamics is governed by a 
fractional Fokker-Planck equation with a spatially varying, periodic potential. As a 
main result we derive a generalization of the celebrated Lifson- Jackson result for normal 
diffusion [HIl ES (HB] to our case with anomalous fractional diffusion: It relates the 
effective fractional diffusion coefficient Ka in Eq. P^ to the bare fractional diffusion 
coefficient in terms of two inverse quadratures of the periodic potential only. As a 
consequence, we find that like in the case with normal diffusion the effective anomalous 
diffusion becomes always suppressed over the bare value. This result may find ample 
application in diverse areas where anomalous diffusion occurs; typical examples are the 
case of superionic conductors [Sn] or for the Josephson junction dynamics PHI ED] when 
the role of disorder may change the normal diffusion into anomalous one. 

In addition, we contrasted the time evolution for normal diffusion with anomalous, 
fractional diffusion. In doing so, we find that after a proper re-scaling of time the 
corresponding asymptotic densities P(x, f) for the coordinate x match each other. 
Distinct differences occur, however, at small evolution times. This time evolution of 
the densities drastically changes upon the application of a finite bias F. Now, the long 
time evolution between normal diffusion and anomalous diffusion becomes markedly 
distinct as well: While the maximum of the biased normal diffusion moves with the 
normal, directed current, the anomalous case is dominated by a ballistic spreading that 
leaves the maximum of the density around the origin. Moreover, this characteristic 
difference can be put to work to differentiate between fractional and fractal Brownian 
diffusion. 
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